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The propagator of a physical degree of freedom ought to obey a Kallen-Lehmann spectral repre- 
sentation, with positive spectral density. The latter quantity is directly related to a cross section 
based on the optical theorem. The spectral density is a crucial ingredient of a quantum field the- 
ory with elementary and bound states, with a direct experimental connection as the masses of 
the excitations reflect themselves into (continuum) 5-singularities. In usual lattice simulational 
approaches to the QCD spectrum the spectral density itself is not accessed. The (bound state) 
masses are extracted from the asymptotic exponential decay of the two-point function. Given the 
importance of the spectral density, each nonperturbative continuum approach to QCD should be 
able to adequately describe it or to take into proper account. In this work, we wish to present a 
first trial in extracting an estimate for the scalar glueball spectral density in SU(3) gluodynamics 
using lattice gauge theory. 
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1. Introduction and Motivation 

The observed hadronic physical spectra can be understood as a collection of bound states of 
quarks and gluons, the fundamental fields of QCD. So far, all attempts to detect directly quarks and 
gluons have produced negative results. However, if one takes QCD seriously, one has to explain 
why quarks and gluons do not contribute to the S-matrix or, from another point of view, why 
quarks and gluons only come about when they are confined into color singlet states. Investigating 
the spectral representation of the QCD propagators does not only help us to perceive better the 
dynamics of the Yang-Mills theory but it also provides clues to grasp the problem of confinement. 

Let &{p) be the Euclidean momentum creation operator associated with a given set of quantum 
numbers. Its vacuum expectation value (&{p) &{—~p)) defines a propagator which encodes all the 
information about the quanta with the chosen quantum numbers. Furthermore, let us assume that 
the propagator has a Kallen-Lehmann spectral representation, i.e. that 

(0(p) ff(-p)) = [ + °°dii , (1.D 

Jo p +H 

where p(ju) is the spectral density, which is real and positive defined if the quantum associated 
with 0{p) contributes physically to the S-matrix. Then, if the spectral density vanishes or becomes 
negative it is a sign that the corresponding quanta are not asymptotic states of the S-matrix. Of 
course, this is not a proof of confinement, but it is a necessary condition to be satisfied by a confined 
particle (assuming that its propagator has a Kallen-Lehmann representation at least). If the quanta 
associated with G contribute to the S-matrix, a typical spectral density is given by 

N 

P (li) = J £z i 8(n-iii) + 0(/x-juo)/(/x), (L2) 

i=\ 

where jJLi are the particle masses squared associated with the operator 0, Z, the corresponding 
probability of creating such a state and \i§ is the threshold where the multiparticle spectrum (= the 
continuum) starts. 

The lattice approach to QCD provides a way to compute nonperturbatively the propagators 
of quarks, gluons and composite operators, such us glueballs, from first principles. The problem 
we address here is what can be learned about the spectral density given a propagator. In particular 
in this contribution, we will revisit the Landau gauge gluon propagator and report on preliminary 
data for the glueball propagator. In the present report we will consider only pure gauge Yang-Mills 
SU(3). For the generation of the gauge configurations, we use the MILC code |jl|]. The simulations 
were performed both at Coimbra and Ghent Universities. 

From the numerical point of view, the computation of p(jU) from the propagator requires an 
inverse integral transform which, typically, calls for some regularization of the data as provided 
by, for example, the maximum entropy method [01. Spectral forms of correlation functions are 
widely studied using lattice simulations, in particular mesons, charmonia etc. at finite T. Note that, 
although it will be interesting to access the propagators in the entire complex p 2 plane, lattice QCD 
simulations only allow us to consider the Euclidean momenta region p 2 > 0. 
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2. On the gluon propagator 

The Landau gauge gluon propagator 

Df v (q) = 8 ab ^ v -^D(q 2 ) (2.1) 

computed on the lattice has been thoroughly scrutinized in the recent years, see for example [||, 
D> H 0, §1- It turns out that the propagator is well described by the tree-level prediction of the 
refined Gribov-Zwanziger action [^] 

rp- _|_ yp- 

"^^TTPWrf <22) 

up to momenta q « 4 GeV with Z ~ 0.81, M 2 ~ 4.0 GeV 2 , M 2 + m 2 ~ 0.57 GeV 2 and A 4 ~ 
0.39 GeV 4 . This functional form means that the low energy propagator has two complex poles 
at 0.412 ± /0.674 GeV. Given that ( |2.2| ) is a perturbative result it is conceivable that taking 
into account quantum corrections will modify the observed behaviour, leading to branch cuts in 
the complex plane. A recent computation of D(p 2 ) over the entire complex p 2 plane using the 



Schwinger-Dyson equations (SDE) [ |10| ] does not observe any complex conjugate poles or other 
cuts than on the negative Euclidean axis. Of course, in order to solve (numerically) the SDE, 
truncations had to be made together with a parametrization of the vertices. All the approximations 
require independent confirmation and further studies are needed to arrive at a reliable conclusion. 



In [11], the gluon spectral function was also probed using fitting forms for the propagator. 

Once the propagator D(q 2 ) is known, one can compute the Schwinger function defined as 



A(0 = J__ ^e-^D{q 2 ) = J dy p(y 2 ) for t > . (2.3) 

If p(p) displays a 8-f unction singularity, then correspondingly A(t) ~ e~ mt and m would be the 
physical mass of an asymptotic free gluon. On the other hand, if the Schwinger function has a 
region where it goes negative, then p(p) must be negative over a given region of p and, therefore, 
the gluon cannot be a free particle in the usual sense. In this latter case we speak about positivity 
violation by the gluon. This violation of positivity has been observed both in lattice calculations 



[ |12[ \13[ |14J and for the solutions of the Schwinger-Dyson equations J10[ , |15J. For completeness, 
in Fig. [I] we show the Schwinger function computed using the lattice gluon propagator for several 
ensembles. Despite the large statistical errors, the plot points towards a p (ju ) with positive and neg- 
ative values and one may conclude that the gluon is, certainly, not an asymptotic state contributing 
to the S -matrix. 



That the gluon is not a physical state can also be viewed from (2.2). However, if one com 



bines the gluon propagator of (Oh to build a glueball propagator [16], at least part of the glueball 



propagator looks like a physical particle propagator and it turns out that its mass value is within the 



values predicted by quenched lattice QCD [17, 18]. In future work, we shall investigate the gluon 



spectral representation using lattice data, based on the hereafter to be explained technology. 
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Figure 1: The Schwinger function A(f ) for the gluon propagator computed for large lattices and several 
volumes and lattice spacings. Details about statistics and lattices sizes and space and be found in m. 



3. On the glueball propagator 



In this section, we will report on an ongoing computation of the spectral density for the color 
singlet glueball operator with quantum number J PC = ++ . In the continuum, such states are 
associated with G = F 2 V . A lattice version of F^ v can be found in Our lattice operator for 
the non-Abelian tensor is such that it reproduces the continuum F^ v up to corrections of order a 2 , 
where a is the lattice spacing. 

The lattice simulations were performed using the Wilson action at j8 = 6.2, meaning a lattice 
spacing a = 0.0726 fm as measured from the string tension, and a lattice of 64 4 with a physical 
size of ~ 4.65 fm. We have generated about 872 gauge configurations. 

For the glueball propagator one expects 



-q)) = infinities polynomial in q 2 + 



-<?)}fi 



mte 



(3-1) 



In order to compute the spectral density, one first has to remove the polynomial infinities. Di- 
mensional analysis shows that the infinities should show up in the coefficients of a polynomial 
^inf(q 2 ) of order q 4 [16]. Therefore, to subtract &j n f(q 2 ) we fit the bare glueball propagator to 
cio + a\q 2 + a\q" over several momenta interval between 4 and 8 GeV and afterwards we subtract 
&i n f(q 2 ) from the bare propagator. For the polynomial fit, in all cases, the reduced % 2 jd.o.f . w 1. 
In Fig. U we report on the bare glueball propagator before subtraction (left). The r.h.s. plot shows 
the subtracted and rescaled glueball propagator defined in such a way that D(q m j n ) = 1, where q m i n 
is the smallest nonvanishing momenta. 
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We would like to call the reader's attention that no data points are shown for q = in Fig. ||. 
The numerical simulations all give 0(0) very large and this point was thence removed from the 
plot. In the continuum, we indeed have that 

D(0) = / + d\i — ► +00 (3.2) 
Jo jU 

since the p(ju)/jU integral also diverges, in good agreement with the discussion of the previous 
paragraph. 




Figure 2: Glueball propagator: (left) bare propagator before subtraction; (right) rescaled subtracted propa- 
gator for various fitting intervals. See text for details. 

The subtracted and rescaled propagator plot show relative independence of the subtraction 
point T, defined as the square of the centre of the fitting interval, indicative that eventual physical 
results will be relatively T-independent. This independence of T allows us to take, from now on, 
only one of the data sets after subtraction. In the following, we will only consider the subtracted 
data obtained from taking the polynomial fit to the range 5 — 6 GeV. 

The finite part, i.e. the glueball propagator after subtraction of ^i n f(q 2 ), is given in terms of 
the spectral density as 

mi) = «%M-„)> flnite = - (i - r) 3 I *"dn 

=-tf-t? r*4i&- 0-3) 

Jo q l + n 

Formally, this equation can be written as follows 

^)=-^=r«„w =J8W)i (3 . 4) 

(q 2 -Ty Jo r+M 

i.e. the l.h.s. of the equation is given by the double Laplace transform of p(ju), where the Laplace 
transform reads ^ fit) = J Q + °° e~ s1 f(s). The propagator is generated by the Monte Carlo technique 
and, therefore, we end up with noisy data as usual. The direct solution of the linear system obtained 
from the last equation is meaningless: it needs to be regularized. Inverting the Laplace transform 
is a typical example of an ill-posed problem. 



5 



Glueball spectral density 



O. Oliveira 



We replaced the original linear system by a Tikhonov-Morozov regularized [19] linear system, 
where the so-called normal equation to be solved is 

J^ 4 p(z) + Ap( z ) = ^(z), 



viz. 



z + t 



(3.5) 



(3.6) 



where A is the regularization parameter. For solving this equation, the (convergent) integrals were 
numerically handled via the introduction of a cutoff at A = 11 GeV and the system was discretized 



via a half-open Newton-Cotes quadrature rule. The solution of ( 3.6) as a function of the momentum 
is show in Fig. ||. 
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Figure 3: p as a function of the momentum. 



The function p(ju) shows a peak around ~ 1.7 GeV and second smaller and broader peak just 
above 3 GeV. These preliminary results compare well the predictions of [17, 18]. There, (contin- 
uum extrapolated) 0++ states with masses of 1710(50)(80) MeV and 2670(180)(130) MeV were 
reported. 
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